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Abstract 

Heat  transfer  characteristics  of  the  metal  hydride  vessel  based  on  the  plate-fin  type  heat  exchanger  were  investigated.  Metal  hydride 
beds  were  filled  with  AB2  type  hydrogen-storage  alloy’s  particles,  Tio.42Zro.58Cro.78Feo.57Nio.2Mno.39Cuo.031  with  a  storage  capacity  of 
0.92  wt.%.  Heat  transfer  model  in  the  metal  hydride  bed  based  on  the  heat  transfer  mechanism  for  packed  bed  proposed  by  Kunii  and 
co-workers  is  presented.  The  time-dependent  hydrogen  absorption/desorption  rate  and  pressure  in  the  metal  hydride  vessel  calculated  by 
the  model  were  compared  with  the  experimental  results.  During  the  hydriding,  calculated  hydrogen  absorption  rates  agreed  with  measured 
ones.  Calculated  thermal  equilibrium  hydrogen  pressures  were  slightly  lower  than  the  measured  hydrogen  pressures  at  the  inlet  of  metal 
hydride  vessel.  Taking  account  of  the  pressure  gradient  between  the  inlet  of  metal  hydride  vessel  and  the  metal  hydride  bed,  it  is  considered 
that  this  discrepancy  is  reasonable.  During  the  dehydriding,  there  were  big  differences  between  the  calculated  hydrogen  desorption  rates  and 
measured  ones.  As  calculated  hydrogen  desorption  rates  were  lower  than  measured  ones,  there  were  big  differences  between  the  calculated 
thermal  equilibrium  hydrogen  pressures  and  the  measured  hydrogen  pressures  at  the  inlet  of  metal  hydride  vessel.  It  is  considered  that 
those  differences  are  due  to  the  differences  of  the  heat  transfer  characteristics  such  as  thermal  conductivity  of  metal  hydride  particles  and 
porosity  between  the  assumed  and  actual  ones.  It  is  important  to  obtain  the  heat  transfer  characteristics  such  as  thermal  conductivity  of 
metal  hydride  particles  and  porosity  both  during  the  hydriding  and  dehydriding  to  design  a  metal  hydride  vessel. 
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1.  Introduction 

Energy  storage  can  play  an  important  role  in  the  opera¬ 
tion  of  electric  energy  systems  ranging  from  load  leveling 
[1]  to  enabling  the  use  of  intermittent  renewable  energy 
sources  such  as  wind  or  solar  energy  for  remote  sites 
to  provide  power  on  demand  [2],  Energy  storage  system 
constructed  around  hydrogen  is  the  combination  of  water 
electrolyzer,  fuel  cell  and  hydrogen- storage  device.  Using 
a  water  electrolyzer,  electricity  production  in  excess  of 
demand  is  converted  to  hydrogen.  Electricity  demand  in 
excess  of  production  is  met  by  converting  hydrogen  to  elec¬ 
tricity  through  a  fuel  cell.  There  are  three  basic  methods 
for  storing  hydrogen:  (1)  compressed  gas;  (2)  liquid;  and 
(3)  chemically  bound  (such  as  metal  hydrides). 

When  using  the  metal  hydrides,  hydrogen- storage  densi¬ 
ties  greater  than  those  achieved  by  liquefaction  are  possible 
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with  gas  pressures  of  just  0. 1  or  0.2  MPa.  The  use  of  metal 
hydrides  also  offers  potential  safety  advantages  over  the 
compressed  hydrogen  or  liquefied  hydrogen.  A  ruptured  ves¬ 
sel  of  compressed  hydrogen  or  liquefied  hydrogen  presents 
an  immediate  and  serious  hazard.  If  the  vessel  ruptured  while 
in  use,  removal  of  the  heat  source  driving  the  thermal  de¬ 
composition  should  cause  the  reaction  to  quench  itself.  In 
addition,  most  metal  hydrides  do  not  react  violently  with  air 
or  water. 

An  a-phase  solid  solution  reacts  with  a  hydrogen  to  form 
(3-phase  hydride  as  shown: 

a-phase  (solid)  +  hydrogen  (gas)  ++  (3-phase  (solid)  (1) 

The  hydriding  reaction  is  normally  exothermic  when  it  is 
proceeding  to  the  right.  On  the  other  hand,  the  dehydrid¬ 
ing  reaction  is  endothermic  when  it  is  proceeding  to  the 
left.  The  performance  of  hydrogen  absorption  to  and  desorp¬ 
tion  from  metal  hydride  vessel  depends  on  the  heat  transfer, 
mass  transfer  and  reaction  kinetics  of  metal  hydride  bed.  Af¬ 
ter  a  few  hydriding  and  dehydriding  cycles,  metal  hydride 
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Nomenclature 

df  pitch  of  fins  (m) 
dp  diameter  of  metal  hydride  particles  (m) 
h*  contact  heat  transfer  coefficient  between  the 
metal  hydride  particles  (W  m-2  K-1) 

/imh  overall  heat  transfer  coefficient  of  the 
metal  hydride  bed  (W  m-2  K-1) 
firs  thermal  radiative  heat  transfer  coefficient 
between  the  metal  hydride  particles 
(W  m-2  K"1) 

firv  thermal  radiative  heat  transfer  coefficient 
between  the  vacant  spaces  (W  m-2  K-1) 
fiw  apparent  heat  transfer  coefficient  in  the 
vicinity  of  bed’s  wall  (W m-2  K-1) 
j  j  factor 

ke  effective  thermal  conductivity  of  metal 
hydride  bed  accompanied  with  the 
hydrogen  flow  (W  m-1  K-1) 
o  effective  thermal  conductivity  of  metal 
hydride  bed  (W  m-1  K-1) 
fiewo  effective  thermal  conductivity  in  the  region 
of  dp/2  from  the  inner  surface  of 
bed’s  wall  (Whi^KT1) 
fif  thermal  conductivity  of  hydrogen 

(W  m-1  K-1) 

fis  thermal  conductivity  of  metal  hydride 

particles  (W  m_1  K_1) 

L\  spacing  between  the  particle  layers  (m) 

L.2  thickness  of  the  vacant  space  (m) 

L3  thickness  of  the  particle  layer  (m) 

Pr  Prandtl  number 

<2i  heat  transfer  rate  through  the  hydrogen 
phase  (W  m-2) 

Q2  heat  transfer  rate  through  the  metal 
hydride  phase-hydrogen  phase-metal 
hydride  phase  (W  m-2) 

Re  Reynolds  number 

AT  temperature  difference  (K) 

Ax  distance  between  the  center  of  solid  layers  (m) 
X  hydrogen  content  (number  of  hydrogen 

atom/number  of  metal  atom) 

Greek  letters 
a  coefficient 

P  coefficient 

y  coefficient 

5  coefficient 

tpi  porosity  in  the  case  where  spheres  are 
packed  loosely 

(p2  porosity  in  the  case  where  spheres  are 
packed  most  compactly 
<py  porosity  of  metal  hydride  bed 
</>w  porosity  of  metal  hydride  bed  in  the 
vicinity  of  bed’s  wall 


particles  break  down  into  fine  particles  of  under  the  several 
lOpm  size.  These  fine  particles  are  good  because  a  large 
surface  area  enhances  both  the  heat  transfer  and  mass  trans¬ 
fer.  However,  fine  particles  leads  to  excessive  pressure  drop 
which  reduces  the  mass  transfer  and  to  small  effective  ther¬ 
mal  conductivity  of  metal  hydride  beds  which  reduces  the 
heat  transfer.  Murthy  and  Gopal  reviewed  the  state-of-the-art 
on  heat  and  mass  transfer  in  the  metal  hydride  beds  [3]. 

Experimental  and  theoretical  studies  of  effective  thermal 
conductivity  in  metal  hydride  beds  have  been  carried  out 
[4—9].  Pons  and  Dantzer  demonstrated  that  small  diameters 
of  metal  hydride  powders  (1 — 40  |xm)  are  responsible  for  the 
low  effective  thermal  conductivity  values  because  of  Knud- 
sen  effect  [10].  They  also  pointed  out  that  the  wall  heat  trans¬ 
fer  coefficient  is  an  important  parameter  for  characterizing 
and  optimizing  the  heat  transfer  efficiency  of  metal  hydride 
bed  [11], 

Different  attempts  to  improve  the  effective  thermal  con¬ 
ductivity  of  metal  hydride  beds  have  been  made:  insertion  of 
aluminum  foam  [12],  integration  of  copper  wire  net  structure 
[13],  packing  in  a  multilayer  waved  sheet  structure  [6],  mi¬ 
croencapsulated  metal  hydride  compact  [14],  and  compact¬ 
ing  metal  hydride  powder  with  an  expanded  graphite  [15]. 

For  effective  design  of  metal  hydride  vessel,  a  mathemati¬ 
cal  model  is  needed  which  effectively  describes  the  physical 
and  chemical  process  in  a  metal  hydride  bed.  A  number  of 
pertinent  and  effective  mathematical  models  for  the  process 
have  been  discussed  [16-23]. 

Metal  hydride  vessels  have  been  designed  and  manufac¬ 
tured  [12,24,25].  The  structure  of  almost  vessels  was  based 
on  the  tubular  type  heat  exchanger.  In  this  paper,  heat  trans¬ 
fer  characteristics  of  the  metal  hydride  vessel  based  on  the 
plate-fin  type  heat  exchanger  was  investigated.  A  model  is 
presented  which  describes  heat  transfer  mechanism  in  the 
metal  hydride  bed  based  on  the  plate-fin  type  heat  exchanger. 
The  time-dependent  hydrogen  transfer  rate  and  pressure  in 
the  metal  hydride  bed  calculated  by  the  model  are  compared 
with  the  experimental  results.  From  the  results  of  those  com¬ 
parison,  the  rationality  and  applicability  of  the  employed 
model  is  discussed  to  design  a  metal  hydride  vessel. 


2.  Experimental  procedure 

2.1.  Hydrogen-storage  alloy 

Hydrogen-storage  alloys  are  broadly  classified  into  three 
types  according  to  the  crystal  structure  and  composition. 
The  first  one  is  AB5  type  hydrogen-storage  alloy  repre¬ 
sented  by  LaNis.  The  second  one  is  AB  type  hydrogen- 
storage  alloy  represented  by  TiFe.  The  third  one  is  AB2 
type  hydrogen- storage  alloy  represented  by  TiCr2,  TiMni  5 
and  one  having  C14  Laves  structure.  The  beds  in  the  metal 
hydride  vessel  were  filled  with  AB2  type  hydrogen-storage 
alloy’s  particles,  Tio.42Zro.58Cro.78Feo.57Nio.2Mno.39Cuo.03, 
with  a  storage  capacity  of  0.92  wt.%.  In  this  hydrogen- storage 
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-0-313.15  K 
(absorption) 
-0-313.15  K 
(desorption) 
-*-353.15  K 
(absorption) 
-*-353.15  K 
(desorption) 


Atom  ratio  (Hydrogen/Metal) 

Fig.  1.  Pressure-composition  temperature  diagram  of  Ti0.42Zr0.58Cr0.7sFe0.57Ni0.2Mn0.39Cu0.03- 


alloy,  Zr  is  substituted  for  a  part  of  Ti,  and  Fe,  Ni,  Mn  and 
Cu  are  substituted  for  a  part  of  Cr  on  the  basis  of  TiCr2. 
Each  substituted  element  has  the  following  effects  [26]: 

(1)  Equilibrium  plateau  pressures  become  variable  by  the 
contents  of  substituted  Zr  and  Fe. 

(2)  Fe  increases  hydrogen- storage  capacity. 

(3)  Fe  makes  the  beginning  point  of  “plateau  region”  more 
distinct. 

(4)  Mn  and  Cu  make  the  gradient  of  “plateau  region” 
smaller. 

(5)  Mn,  Ni  and  Cu  make  the  hysteresis  between  the  hydrid- 
ing  process  and  dehydriding  process  smaller. 

A  pressure-composition  temperature  diagram  of  this 
hydrogen- storage  alloy  is  shown  in  Fig.  1.  Absorption  and 
desorption  pressure  plateau  of  the  examined  metal  hydride 
are  0.37  and  0.31  MPa  at  the  hydrogen/metal  ratio  of  0.5 
and  the  temperature  of  313.15  K,  respectively. 

2.2.  Metal  hydride  vessel 

The  hydriding  reaction  in  the  metal  hydride  is  exothermic. 
On  the  other  hand,  the  dehydriding  reaction  is  endothermic. 
Under  hydriding/dehydriding  reactions,  the  metal  hydride 
beds  must  be  heated  or  cooled.  For  this  reason,  the  structure 
of  the  metal  hydride  vessel  is  the  same  as  that  of  the  heat 
exchanger. 

In  general,  owing  to  the  fine  particles  of  metal  hydride,  the 
effective  thermal  conductivity  of  metal  hydride  bed  usually 
has  a  value  around  1  Wm_1  K-1  at  the  ambient  tempera¬ 
ture.  Thus,  the  heat  transfer  resistance  of  metal  hydride  bed 
determines  the  overall  heat  transfer  resistance  in  the  case 
where  liquid  is  used  to  heat  and  cool  the  metal  hydride  beds. 

The  methods  suggested  for  heat  and  mass  transfer  en¬ 
hancement  can  be  broadly  divided  into  two  main  groups.  In 
the  first  method,  the  heat  transfer  effectiveness  of  the  metal 
hydride  bed  is  improved  by  providing  extended  area  for  heat 
transfer  in  the  forms  of  fins,  foams  or  meshes.  Metal  hydride 


particles  are  distributed  into  the  matrix  formed  by  the  high 
thermal  conductivity  materials  such  as  copper  or  aluminum. 
In  the  second  method  of  enhancement,  the  metal  hydride 
particles  are  consolidated  through  a  binder  or  metal  hydride 
particle.  High  thermal  conductivity  materials  such  as  copper 
or  carbon  are  compressed  under  high  pressure  into  porous 
metal  matrix  hydride  compacts. 

The  first  method  by  fins  was  applied  in  this  research.  An 
exterior  view  of  the  metal  hydride  vessel  is  shown  in  Fig.  2. 
The  cross-sectional  view  is  shown  in  Fig.  3.  Specifications 
are  shown  in  Table  1 .  The  concept  of  brazed  plate- fin  type 
heat  exchanger  was  applied  to  extend  the  heat  transfer  area 
between  the  metal  hydride  particles  and  water.  The  vessel 
was  composed  of  an  aluminum  alloy  rectangular  body  and 
includes  three  packed  beds  for  the  metal  hydride  particles 
and  four  passages  for  the  water.  The  metal  hydride  particles 
were  filled  between  the  fins. 

2.3.  Experimental  setup 

The  measurements  for  the  heat  transfer  characteristics  of 
metal  hydride  vessel  were  performed  using  the  experimen¬ 
tal  setup  shown  schematically  in  Fig.  4.  Metal  hydride  ves¬ 
sel  was  provided  with  hydrogen  from  a  high-pressure  cylin¬ 
der  sold  at  a  market.  Hydrogen  in  a  high-pressure  cylinder 
was  a  byproduct  from  sodium  chloride  electrolysis,  which 
was  compressed,  purified  and  filled  up.  Purity  and  dew  point 
of  hydrogen  was  99.998%  and  213.15  K,  respectively.  Im¬ 
purity  was  nitrogen  of  which  concentration  was  19.3  ppm 
(0.00193  vol.%). 

After  maintaining  the  metal  hydride  vessel  at  the  tem¬ 
perature  and  pressure  of  experimental  conditions,  hydrogen 
was  absorbed  into  the  metal  hydride  vessel  or  desorbed  from 
the  metal  hydride  vessel.  The  hydrogen  flow  rate  was  kept 
constant  at  2.1  or  3.3  Nm3  h_1.  Water  of  0.6  or  1.2m3  h-1 
was  delivered  to  the  vessel  to  remove  or  provide  the  reaction 
heat.  The  temperature  distribution  in  the  metal  hydride  bed, 
water  inlet  and  outlet  temperature,  hydrogen  pressure  at  the 
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Fig.  2.  Exterior  view  of  the  metal  hydride  vessel. 


Table  1 

Specifications  of  the  metal  hydride  vessel 


Material 

Hydrogen-storage  alloy 

Tio.42Zro.58Cro.7gFeo.57Nio.2Mno.39Cuo.03 

Vessel 

Aluminum  alloy 

Storage  capacity  of 

>1 

hydrogen  (Nm3) 

Weight  (kg) 

Hydrogen-storage  alloy 

11.3 

Vessel 

19.7 

Total 

31 

Overall  size  (mm) 

Height 

105 

Length 

210 

Width 

326 

Heat  transfer  area  (m2) 

Metal  hydride  beds 

2.85 

Flow  passages  for  water 

0.81 

Fig.  4.  Experimental  setup. 
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inlet  of  metal  hydride  vessel,  hydrogen  flow  rate  and  wa¬ 
ter  flow  rate  were  gathered  using  data  acquisition  equip¬ 
ment. 


3.  Heat  transfer  mechanism  in  the  metal 
hydride  bed 

It  can  be  considered  that  heat  transfer  mechanism  in  the 
metal  hydride  bed  is  the  same  as  that  in  the  common  pow¬ 
der  packed  bed  such  as  catalyst  bed.  Kunii  and  co-workers 
discussed  the  heat  transfer  mechanism  in  packed  bed  using 
the  effective  thermal  conductivity  in  packed  bed  and  appar¬ 
ent  heat  transfer  coefficient  in  the  vicinity  of  vessel’s  wall 
[27-30].  With  reference  to  those  discussions  heat  trans¬ 
fer  mechanism  in  the  metal  hydride  bed  can  be  described 
as  follows,  provided  that  the  following  assumptions  are 
made: 

(1)  The  metal  hydride  particles  and  hydrogen  in  the  vacancy 
of  metal  hydride  bed  are  in  thermal  equilibrium. 

(2)  There  is  no  pressure  gradient  within  the  metal  hydride 
bed. 

(3)  The  temperature  in  the  metal  hydride  bed  is  only  a  func¬ 
tion  of  time. 

3.1.  Effective  thermal  conductivity  in  the  metal 
hydride  bed 

Yagi  and  Kunii  gave  an  effective  thermal  conductivity  in 
the  packed  bed  accompanied  with  the  fluid  flow  as  follows 
[27]: 

h  =  k^+aPrRe  (2) 

kf  kf 

As  reactions  in  the  metal  hydride  bed  is  usually  fast,  it  can 
be  considered  that  there  is  not  hydrogen  flow  in  the  metal 
hydride  bed.  Even  though  the  whole  quantity  of  hydrogen 
introduced  to  or  discharged  from  the  vessel  does  not  re¬ 


act,  Reynolds  number  based  on  the  diameter  of  metal  hy¬ 
dride  particles  is  very  small.  Thus,  it  can  be  considered  that 
metal  hydride  bed  is  not  accompanied  with  the  hydrogen 
flow. 

Heats  are  transferred  in  parallel  by  the  following  mecha¬ 
nisms  in  the  case  where  temperature  gradient  exists  in  the 
metal  hydride  bed  as  shown  in  Fig.  5: 

(1)  Contact  heat  transfer  between  the  metal  hydride  parti¬ 
cles. 

(2)  Heat  conduction  through  the  thin  hydrogen  film  in  the 
vicinity  of  the  contact  point  of  metal  hydride  particles. 

(3)  Heat  radiation  from  the  surface  to  other  surfaces  of  the 
metal  hydride  particles. 

(4)  Heat  conduction  through  the  metal  hydride  particles. 

(5)  Heat  conduction  through  the  hydrogen  in  a  vacant  space. 

(6)  Heat  radiation  from  the  vacant  space  to  the  adjacent 
another  vacant  space. 

Heat  transfer  mechanism  is  modeled  as  shown  in  Fig.  6. 
Heat  transfer  rate  through  the  hydrogen  phase  per  unit  area 
is  described  as  follows: 

e,  =  (£+*„)  AT  (3) 

Heat  transfer  rate  through  the  metal  hydride  phase-hydrogen 
phase-metal  hydride  phase  per  unit  area  is  described  as  fol¬ 
lows: 

AT 

Ql  =  (1  /{h*  +  (kf/L2)  +  hrs))  +  (L3/ks)  '4S 

It  can  be  considered  that  heat  between  the  metal  hydride 
particles,  vacant  spaces  and  metal  hydride  particles  is  trans¬ 
mitted  through  the  cross-sectional  area  of  metal  hydride  bed 
equivalent  to  1  —  0V.  On  the  other  hand,  it  can  be  considered 
that  that  between  the  vacant  spaces  is  transferred  through 
the  cross-sectional  area  of  metal  hydride  bed  equivalent  to 
</>v  Temperature  gradient  between  the  surfaces  C  and  D  in 
Fig.  6  is  AT/L] .  Effective  thermal  conductivity  that  can  de¬ 
scribe  the  total  heat  transfer  in  the  metal  hydride  bed  is 


Fig.  5.  Heat  transfer  in  the  metal  hydride  bed. 
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introduced  as  follows: 

tt0I7=*’(lT  +  ,"v)A7' 

+ _ _  (3) 

(1  /(h*  +  (M/Ll)  +  Ars))  +  (**/*■) 

L\  can  be  considered  as  the  average  distance  between  the 
centers  of  adjacent  two  metal  hydride  particles  and  is  de¬ 
scribed  as  follows: 


U  =  (6) 

L.2  and  L3  are  described,  respectively,  as  follows: 

Li  =  ydp  (7) 

L3  =  Sdp  (8) 


Substituting  Eqs.  (6)— (8)  into  Eq.  (5),  effective  thermal  con¬ 
ductivity  of  the  metal  hydride  bed  is  described  as  follows: 


+ _ _ 

(1  /my)  +  Wdp/ktHh*  +  *«)))  +  (8/pm/ks) 


In  the  case  where  fluid  exists  in  the  vacant  spaces  of  packed 
bed,  contact  heat  transfer  resistance  between  the  particles  is 
smaller  than  other  heat  transfer  resistances  and  is  negligible. 
When  it  is  assumed  that  metal  hydride  particles  are  spheres 
of  diameter  dp  and  particle  layers  in  Fig.  6  are  rectangular 
bodies  of  which  cross-sectional  area  is  7tdp4  and  volume 
is  jrt/p/6,  S  is  determined  as  2/3.  As  the  metal  hydride  bed 
investigated  in  this  paper  is  operated  around  the  ambient 
temperature,  heat  radiation  is  so  small  that  this  heat  transfer 
mechanism  is  negligible,  ft  is  considered  as  1.  Therefore, 
effective  thermal  conductivity  in  the  metal  hydride  bed  is 


simplified  as  follows: 

AeO  =  l~0v 

kf  ^  y+(2/3)(k{/ks) 


(10) 


y  gives  the  heat  transfer  rate  through  the  fluid  phase  near  the 
contact  surface  of  solid  particles.  Kunii  and  Smith  obtained 
y  of  two  cases  [28].  First  one  is  y\  for  the  case  where  spheres 
are  packed  loosely  and  porosity  is  0.476.  Second  one  is  yi 
for  the  case  where  spheres  are  packed  most  compactly  and 
porosity  is  0.26.  It  is  supposed  that  actual  packed  bed  may 
be  a  mixture  of  these  two  packing  cases  and  y  in  actual 
packed  bed  can  be  described  as  follows: 


y  =  yz  +  (yi  —  yi) 


<ftv  -  <p2 

<Pl  -  (p2 


=  yi  +  (yt  -  yi) 


<pv  -  0.26 
0.216 


(11) 


3.2.  Apparent  heat  transfer  coefficient  in  the  vicinity  of 
bed’s  wall 

Yagi  and  Kunii  introduced  an  apparent  heat  transfer  co¬ 
efficient  in  the  vicinity  of  tube’s  wall  in  the  case  where  a 
fluid  of  high  temperature  is  flowed  into  the  packed  bed  and 
is  cooled  steadily  by  the  tube’s  wall  of  constant  tempera¬ 
ture  [29].  Shape  of  particles  applied  in  their  research  was 
sphere.  It  was  assumed  that  the  ratio  between  the  diame¬ 
ter  of  particles  and  inner  diameter  of  tube  was  sufficiently 
small  and  the  temperature  difference  between  the  fluid  and 
particles  was  too  small  to  be  negligible.  In  this  paper,  those 
assumptions  are  applied. 

It  is  assumed  that  temperature  distribution  in  the  metal 
hydride  bed  is  determined  by  the  effective  thermal  conduc¬ 
tivity  of  metal  hydride  bed  as  described  in  Eq.  (10).  Using 
the  effective  thermal  conductivity  in  the  region  of  dp! 2  from 
the  inner  surface  of  bed’s  wall,  apparent  heat  transfer  coef¬ 
ficient  in  the  vicinity  of  bed’s  wall  is  determined  as  follows: 


kf  {kf/  ke wo)  -  {kf  /  ke 0) 
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In  the  same  manner  as  Eq.  (9),  effective  thermal  conductivity 
in  the  region  of  dp/2  from  the  inner  surface  of  bed’s  wall  is 
described  as  follows: 


+ _ ^ _ 

(1/((1/Xw)  +  C hIsdp/kf )))  +  (kf/3ks) 


(13) 


As  mentioned  earlier,  heat  radiation  is  so  small  that  this 
heat  transfer  mechanism  is  negligible.  Therefore,  Eq.  (13) 
is  simplified  as  follows: 


ke wO  _  ~  ,  1  —  0w 

kf  Kw  +  (kf/3ks) 


(14) 


yw  is  described  as  follows: 

((«  -  \)/k)2  1 

/w  4(lnA:  —  ((/c  —  l)//c))  3k 

_  ks 
K~kf 


(15) 

(16) 


3.3.  Overall  heat  transfer  coefficient  of  the  metal 
hydride  bed 

Heat  transfer  resistance  of  the  metal  hydride  bed  is  the 
sum  of  that  of  central  part  in  the  metal  hydride  bed  and  that 
in  the  vicinity  of  the  bed’s  wall.  Thus,  overall  heat  transfer 
coefficient  of  the  metal  hydride  bed  is  described  as  follows: 

jr—£-  +  r  (it) 

hua  2  keo  hw 


4.  Experimental  results 

Measurement  results  of  time-dependent  temperature  of 
the  metal  hydride  bed  and  pressure  at  the  inlet  of  the 
metal  hydride  vessel  during  the  hydriding  at  the  inlet  water 
temperature  of  293.15  K  and  at  the  hydrogen  flow  rate  of 
3.3Nm3h_1  are  shown  in  Fig.  7.  Those  during  the  dehy- 
driding  at  the  inlet  water  temperature  of  313.15  K  and  at 
the  hydrogen  flow  rate  of  3.3  N  m3  h-1  are  shown  in  Fig.  8. 
The  difference  of  water  flow  rate  does  not  influence  on  the 


0  5  10  15  20 

Time  (min) 


Water  flow  rate:  0.6  mV1 

(1)  Temperature  of  the  metal  hydride  bed 

(2)  Pressure  at  the  inlet  of  metal  hydride  vessel 


Water  flow  rate:  1.2  mV1 

( 1 )  Temperature  of  the  metal  hydride  bed 

(2)  Pressure  at  the  inlet  of  metal  hydride  vessel 


Fig.  7.  Measured  results  during  the  hydriding.  Hydrogen  flow  rate,  3.3Nm3h  *;  inlet  water  temperature,  293.15 K. 


Fig.  8.  Measured  results  during  the  dehydriding.  Hydrogen  flow 


e,  3.3  Nm3  h_1; 


inlet  temperature,  313.15  K. 
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metal  hydride  temperature  and  pressure  at  the  inlet  of  the 
metal  hydride  vessel  during  both  the  hydriding  and  dehy- 
driding.  The  effective  thermal  conductivity  of  metal  hydride 
bed  usually  has  a  value  around  1 W  m-1  K-1  at  the  ambient 
temperature.  Thus,  the  difference  of  water  flow  rate  does 
not  influence  on  the  overall  heat  transfer  resistance  of  metal 
hydride  bed. 

5.  Discussion 

5.1.  Dependences  of  effective  thermal  conductivity 
of  metal  hydride  bed  on  thermal  conductivity  of 
metal  hydride  particles 

Using  these  equations,  calculated  effective  thermal  con¬ 
ductivities  of  metal  hydride  bed  at  293.15  K  are  shown  in 
Fig.  9.  Effective  thermal  conductivities  of  metal  hydride 
bed  increase  with  increase  of  thermal  conductivity  of  metal 
hydride  particles  and  with  decrease  of  porosity  of  metal 
hydride  bed.  Generally,  in  the  actual  bed  of  uniform  and 
small  particles,  porosity  is  in  the  range  between  0.36  and 
0.4  [30].  Measured  effective  thermal  conductivity  of  metal 
hydride  bed  at  the  ambient  temperature  has  a  value  around 
1 W  m-1  K-1  [4,9].  When  it  is  assumed  that  porosity  of  the 
metal  hydride  bed  is  0.4  and  thermal  conductivity  of  metal 
hydride  particles  is  in  the  range  from  5  to  10Wm_1  K-1, 
calculated  effective  thermal  conductivities  of  metal  hydride 
bed  are  almost  equal  to  the  measured  results. 

5.2.  Comparison  between  the  measured  and  calculated 
results 

Calculated  results  of  hydrogen  flow  rate  absorbed  to  or 
desorbed  from  the  metal  hydride  vessel  and  thermal  equilib¬ 
rium  pressure  are  compared  with  the  measurement  results 
in  order  to  verify  the  propriety  of  heat  transfer  mechanism 
in  the  metal  hydride  bed  described  in  the  earlier  discussion 
as  follows. 

It  is  reported  that  the  effective  thermal  conductivity  of 
the  metal  hydride  bed  increases  with  the  absorption  of  hy¬ 


drogen.  The  large  increase  in  effective  thermal  conductiv¬ 
ity  or  decrease,  respectively,  in  the  plateau  region  is  due  to 
the  large  change  of  volume  during  absorption  or  desorption 
of  hydrogen.  So  it  is  assumed  that  thermal  conductivity  of 
metal  hydride  particles  is  varied  with  the  hydrogen  content. 
Alloys  based  on  titanium  have  thermal  conductivity  in  the 
range  of  7-15  Wm_1  K-1  at  the  ambient  temperature  [31]. 
With  reference  to  those  values,  it  was  assumed  that  thermal 
conductivity  of  metal  hydride  particles  was  given  by  the  fol¬ 
lowing  equation: 

ks  =  3  +  5X  (18) 

It  was  assumed  that  porosity  of  metal  hydride  bed  was  0.4 
[30]  and  that  in  the  vicinity  of  bed’s  wall  was  0.7  [29]. 
Heat  transfer  coefficient  of  water  was  calculated  from  the 
following  equation: 

j  =  0.72 Re-°J  (19) 

Hydrogen  content  in  the  metal  hydride  vessel  at  the  be¬ 
ginning  of  hydrogen  absorption  or  desorption  was  deter¬ 
mined  using  the  measured  temperature  and  pressure  and 
the  pressure-composition  temperature  diagram  of  metal  hy¬ 
dride.  Considering  the  enhancement  of  heat  transfer  by  ex¬ 
tended  area  of  aluminum  fins,  temperature  of  vessel’s  wall, 
overall  heat  transfer  coefficient  of  the  metal  hydride  bed  and 
heat  transfer  coefficient  of  water  were  determined  so  that 
heat  transfer  rate  of  metal  hydride  bed  was  equal  to  that  of 
water.  Heat  transfer  rate  divided  by  the  heat  of  hydriding 
reaction  or  that  of  dehydriding  reaction  gave  hydrogen  flow 
rate  absorbed  to  or  desorbed  from  the  metal  hydride  vessel. 
Addition  of  increment  of  hydrogen  content  to  or  subtraction 
of  reduction  of  hydrogen  content  from  that  of  the  previous 
time  gave  the  hydrogen  content  at  that  time.  Thermal  equi¬ 
librium  pressure  of  metal  hydride  bed  at  that  time  was  deter¬ 
mined  using  the  measured  temperature,  obtained  hydrogen 
content  and  pressure-composition  temperature  diagram  of 
metal  hydride. 

Comparisons  between  the  measured  results  and  calculated 
ones  during  the  hydriding  and  during  the  dehydriding  are 
shown  in  Figs.  10  and  11,  respectively.  During  the  hydriding, 


Thermal  conductivity  of  metal  hydride  particles  (Wm't  ) 


Fig.  9.  Dependences  of  effective  thermal  conductivity  of  metal  hydride  bed  on  thermal  conductivity  of  metal  hydride  particles. 


313.15  K. 


calculated  hydrogen  absorption  rates  agreed  with  measured 
ones.  Calculated  thermal  equilibrium  hydrogen  pressures 
were  slightly  lower  than  the  measured  hydrogen  pressures 
at  the  inlet  of  metal  hydride  vessel.  Taking  account  of  the 
pressure  gradient  between  the  inlet  of  metal  hydride  vessel 
and  the  metal  hydride  bed,  it  is  considered  that  those  dis¬ 
crepancies  are  reasonable. 

In  this  paper,  it  is  assumed  that  there  are  no  differences 
of  the  heat  transfer  characteristics  such  as  thermal  con¬ 
ductivity  of  metal  hydride  particles  and  porosity  between 
during  the  hydriding  and  during  the  dehydriding.  However, 
for  example,  the  effective  thermal  conductivity  during  the 
dehydriding  is  higher  than  that  during  the  hydriding  [9]. 
This  hysteresis  is  probably  due  to  the  fact  that  the  hydrogen 
content  during  the  dehydriding  is  somewhat  higher  than 
during  the  hydriding  at  the  same  pressure.  As  a  result,  dur¬ 
ing  the  dehydriding  there  were  big  differences  between  the 
calculated  hydrogen  desorption  rates  and  measured  ones. 
As  calculated  hydrogen  desorption  rates  were  lower  than 
measured  ones,  there  were  big  differences  between  the  cal¬ 


culated  thermal  equilibrium  hydrogen  pressure  and  the  mea¬ 
sured  hydrogen  pressure  at  the  inlet  of  metal  hydride  vessel. 


6.  Conclusions 

Heat  transfer  characteristics  of  the  metal  hydride  vessel 
based  on  the  plate-fin  type  heat  exchanger  were  investigated. 
Calculated  results  using  the  heat  transfer  model  in  the  metal 
hydride  bed  based  on  the  heat  transfer  mechanism  for  packed 
bed  proposed  by  Kunii  and  co-workers  were  compared  with 
the  measured  results.  During  the  hydriding,  calculated  hy¬ 
drogen  absorption  rates  agreed  with  measured  ones.  Calcu¬ 
lated  thermal  equilibrium  hydrogen  pressures  were  slightly 
lower  than  the  measured  ones  at  the  inlet  of  metal  hydride 
vessel.  Taking  account  of  the  pressure  gradient  between  the 
inlet  of  metal  hydride  vessel  and  the  metal  hydride  bed,  it 
is  considered  that  those  discrepancies  are  reasonable. 

There  are  differences  of  the  heat  transfer  characteristics 
such  as  thermal  conductivity  of  metal  hydride  particles  and 
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porosity  between  during  the  hydriding  and  during  the  de- 
hydriding.  As  a  result  during  the  dehydriding  there  were 
big  differences  between  the  calculated  hydrogen  desorption 
rates  and  measured  ones.  As  calculated  hydrogen  desorp¬ 
tion  rates  were  lower  than  measured  ones,  there  were  big 
differences  between  the  calculated  thermal  equilibrium  hy¬ 
drogen  pressures  and  the  measured  ones  at  the  inlet  of  metal 
hydride  vessel. 

It  is  important  to  obtain  the  heat  transfer  characteristics 
such  as  thermal  conductivity  of  metal  hydride  particles  and 
porosity  both  during  the  hydriding  and  dehydriding  to  design 
a  metal  hydride  vessel. 
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